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We re- visit the phase diagram of the frustrated spin-1 /2 ladder with two competing inter-chain antiferromag- 
netic exchanges, rung coupling J± and diagonal coupling J x . We suggest, based on the accurate renormal- 
ization group analysis of the low-energy Hamiltonian of the ladder, that marginal inter-chain current-current 
interaction plays central role in destabilizing previously predicted intermediate columnar dimer phase in the 
vicinity of classical degeneracy line Jj_ = 2 J x . Following this insight we then suggest that changing these 
competing inter-chain exchanges from the previously considered antiferromagnetic to the ferromagnetic ones 
eliminates the issue of the marginal interactions altogether and dramatically expands the region of stability of 
the columnar dimer phase. This analytical prediction is convincingly confirmed by the numerical density matrix 
renormalization group and exact diagonalization calculations as well as by the perturbative calculation in the 
strong rung-coupling limit. The phase diagram for ferromagnetic J± and J x is determined. 

PACS numbers: 75.10Jm, 75.10.Pq, 75.30.Kz, 75.40.Cx 



I. INTRODUCTION 

Frustrated quantum antiferromagnets have for a long time 
attracted attention of both theorists and experimentalists^ One 
of the main reasons for this continuing focus is the dominant 
role of quantum fluctuations in stabilizing various classical 
and quantum orders in this class of systems. Spin ladders 
represent particularly interesting class of models exhibiting 
rich variety of phases 2 - - — . In addition to being interesting in 
their own right, spin ladders allow for high precision numer- 
ical investigations by density matrix renormalization group 
(DMRG) 7,8 and quantum Monte Carlo techniques. One of the 
most striking features of spin ladders consists in the finding 2 
that generic inter-chain interaction flows toward strong cou- 
pling, resulting in a confinement of gapless spin- 1/2 spinon 
excitations of constituent spin chains (which form legs of the 
ladder) into tightly bound spin-1 pairs (triplons). Extensive 
experimental efforts 9 resulted in observation of one- and even 
two-triplons states in La4Si"ioCu2404i 10 . Most recently, an 
evolution of spin excitations from those of deconfined spinons 
at high energies to bound triplet and singlet spinon pairs at 
low energies has been mapped via inelastic neutron scattering 
in the weakly coupled ladder material CaC^Ojii. 

While the standard ladder geometry realizes inter-chain 
interaction between spins on the legs in the form of non- 
frustrated exchange along the rungs of the ladder, a more com- 
plicated geometry is possible as well. In this work we focus 
on the ladder with frustrated inter-chain interactions between 
leg spins, when the inter-chain exchange takes place both 
on rungs (equation ([3]) below) and diagonals (equation (O). 
Such geometry is in fact typical for many spin chain materi- 
als where the superexchange between spins proceeds via 90- 
degree Cu-O bonds. In particular, well studied spin chain ox- 
ides SrCu02 and S^CuOs are characterized by the presence 
of (very weak) rung and diagonal exchanges 12 which frustrate 
correlations between chain spins and leads to extremely low 
three-dimensional ordering temperatures. 

In this work, we re-visit and resolve one of the outstand- 
ing questions in this field - the appearance of spontaneously 



dimerized ground state in frustrated spin ladder model with 
only pairwise exchange interactions between microscopic lat- 
tice spins. Specifically, the Hamiltonian of the problem reads 



H = H U 



H, 



diag 



where 



H\ cg — J fSl,n ' Si, n +i + S2,n ' S 



2,n+l 



(1) 



(2) 



describes two isotropic Heisenberg chains with positive (anti- 
ferromagnetic) nearest-neighbor exchange J while 



and 



diag 



-Hrung — J± Sl,n 1 S2,; 



(3) 



(4) 



describe frustrated inter-chain interaction H- mtcr = ^rung + 

-*^diag- 

In the weak-coupling limit {J±,J X <C J) one treats inter- 
chain interaction H mtcl as a perturbation and takes contin- 
uum limit along the chain direction. It is then easy to see 
that for J± < 2 J x the ground state is of Haldane type, with 
two spin- 1/2 on the rung forming effective spin-1, while for 
J± > 2J X rung pairs form singlets, resulting in the rung- 
singlet (RS) phase. Transition region J± m 2 J x between 
these two well-understood phases requires careful analysis 
which is described in Reffl3l It was found there that in the 
narrow region (boundaries are approximate) 



2J X 
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< J± < 2 Jx 
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the ladder should exhibit columnar dimer (CD) phase in its 
ground state 13 . This finding was questioned in several exten- 
sive numerical studies 1415 which suggested that the CD phase 
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is absent and that instead there is a direct transition between 
Haldane (J± < 2J X ) and RS (J± > 2J X ) ground states. 
The most recent work~ on this subject does find the dimer- 
ized phase, although the evidence for this is not particularly 
strong. 

Narrow extent (O of the suggested CD order makes numer- 
ical analysis of the problem difficult. We will argue below 
that in the case of antiferromagnetic couplings, J±,J X > 0, 
the situation is even more complicated by the presence of 
marginally relevant inter-chain interaction between spin cur- 
rents (uniform magnetization) of the two chains. We show that 
this interaction is responsible for suppressing the CD instabil- 
ity for not too small inter-chain exchange values (J± > 0.3 J 
or so) and producing the first order phase transition between 
the Haldane and RS phases instead. We also show that chang- 
ing the sign of the inter-chain couplings to a ferromagnetic 
one (so that Jj_ , J x < 0) effectively removes the current- 
current interaction from the problem and allows one to access 
the CD phase even for not too small |Ji. x | values. These 
arguments, derived from renormalization group (RG) analy- 
sis described in Section|II] are supported by extensive DMRG 
and exact-diagonalization calculations as well as the pertur- 
bation analysis for the strong rung-coupling limit, results of 
which are reported in Section [Til] The ground-state phase di- 
agram in the J± - J x plane for the ferromagnetic case is also 
presented there. The case of antiferromagnetic couplings be- 
tween chains of the ladder is addressed again in Section [IV] 
We conclude by summarizing our findings in Section [V] Ap- 
pendix describes an application of our numerical approach 
to the well-understood case of a single frustrated Heisenberg 
chain which is known to realize the spontaneously dimerized 
ground state. 



II. RG ANALYSIS 

Low-energy description of the problem is based on the con- 
tinuum representation of the spin operator 



S n /a =M(x) + (-l) x N(x), 



(6) 



where x = nag and clq is the lattice spacing, which we 
set to unity in what follows. Uniform M and staggered N 
magnetizations represent spin fluctuations with momenta near 
q = and ir/ao, correspondingly. Another very important 
for the following low-energy degree of freedom is staggered 
dimerization e(x). It represents fluctuational part of the bond 
strength (energy density) between two neighboring spins on 
the chain, 



e(x) = (-l) x (s n -S nH 



-a 



(7) 



The second term, G av , represents an average energy per bond, 
which is a position-independent constant. Spontaneously 
dimerized ground state is characterized by the finite expecta- 
tion value of dimerization, (e), which describes the staggered 
pattern of strong and weak bonds along the chain. 



Low-energy limit of the interchain Hamiltonian is then 
found to contain at least 5 couplings that flow under RG. 

Winter = / dx (j)iN + g 2 J + g 3 a 2 M + g±E + g 5 K^j (8) 
where we introduced the following short-hand notations 



JV = Ni -N 2 , J = M 1R -M 2L 
M = d x N ± ■ d x N 2 , E = eie 2 , 
K = Mijj • Mil + M 2R ■ M 2L 



M 1L -M 2R , 



(9) 



The relevant couplings are 171,4 which describe coupling be- 
tween staggered magnetizations Nj and staggered dimeriza- 
tions ej of the chains, correspondingly. Marginal couplings in- 
clude g 2 and g$ which describe current-current interaction be- 
tween chains (g 2 ) as well as residual (and naively, marginally 
irrelevant) in-chain backscattering (55). The irrelevant terms 
contain 93, which is of key importance for generating (to- 
gether with g 2 ) the novel £162 term (coupling 94 above). 
This term should be, strictly speaking, be supplemented with 
ged x eid x e 2 which will certainly be generated (as well as other 
more irrelevant terms) - but its small initial value, ~ J\ x /«/, 
allows us to neglect it. Initial values of the couplings are: 

.91 (0) = J± - 2J X , 0,(0) = J± + 2J X , 53 (0) = Ji/2, 
54 (0) = , g 5 (0) = -0.23(2™). (10) 

The value of g$ (0) has been estimated in Ref. [T^. RG equa- 
tions for our mode l 13 ' 15 are rather similar to the much studied 
case of the zig-zag ladder 18 . They are conveniently formu- 
lated in terms of dimensionless variables 

Gl , 3A = i^,G 2 , 5 = p!> (11) 
Zirv Zirv 

and read (g = dg/dl, where I is the RG scale) 

G\ = G± + G\G 2 — -G1G5 + G2G3 — -G2G4, 

G 2 = G^+4G?+4GiG 3 -4GiG4-2G 3 G4, 
G3 = — G3, 

3 3 3 

G4 = G4 — -G1G2 — -G2G3 + -G4G5, 



G5 — G^ — 2G^ + 2G1G3 2G4 



(12) 



Almost everywhere in parameter space (J±, J x ) 31,4 di- 
verge exponentially while 173 dies off exponentially and the 
flow is completely determined by the initial value of <?i(0) 
(remember that 174(0) = 0). But in the vicinity of J± = 2 J x 
line things are different because naively there <?i(0) = as 
well. Thus there two relevant operators appear to be absent 
due to fine-tuning. This region requires careful study. 

Let us assume that we are very close to this line so that 
Gi(0) = 0(J x /J 2 ). Note that other couplings are G 2 , 3(0) = 
0(J x /J) because they are not sensitive to the difference 
J± — 2 J x . In order to understand how the flow starts consider 
very short RG times, £ < £q = 1 and keep only terms of order 
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(Jx/J) 2 in the right-hand-side of the RG equations. The sys- 
tem simplifies dramatically and we find that G 2! 3.5 'decouple' 
from the equations for the relevant G14 couplings 

G 2 (l) = G 2 (0)/(1-G 2 (0)f), 
G 3 (£) = G 3 (0)e- e , 

G 5 (£) = G B (0)/(l + |G B (0)|*)->-l/*. (13) 
Thus for £ ~ 1 we have 

G\ = Gi + G2G3, 

G 4 = G 4 -^G 2 G 3 . (14) 

Note that for such short £'s we can safely treat marginal cou- 
plings G 2 ,5 as constants, because G 2 .s(0) <C 1. The last two 
equations then admit straightforward solution 

GiW = (Gi(0) + fy-fe-*, 

Gi(t) = ^sinh^], (15) 

where q = G 2 (0)G 3 (0) = {2J x /2nv) 2 = (2J x /tt 2 J) 2 , and 
we made use of the initial condition G4(0) = and of the 
fact that the spin velocity is given by v = ttJ/2. We ob- 
serve renormalization of the initial values of G\^ couplings 
by quantum fluctuations. It is clear that for intermediate range 
of < £ <C £\ = ln( J I J x ), where the approximation < TT~4T > 
is expected to work, the two equations ( fT51 ) describe evolu- 
tion of relevant G\^ couplings with initial values given by 
Gi(0) -> Gi(0) + q/2 and G 4 (0) ->■ -3<?/4, correspond- 
ingly. The same renormalization can be obtained via direct 
perturbative calculation, as was done originally in Ref.ll3l 

Once Gi ; 4 got some non-zero initial values of the order 
(J x /J) 2 , they will grow exponentially fast and at £\ ~ 
ln(J/J x ) < £ 2 = 1/G 2 (0) — J/J x reach values of the 
order J x /J. Note that G 2 (ii) ~ Jx/J < 1 at this scale 
and our weak-coupling consideration still makes sense. An 
estimate of the range where the dimerized phase is expected 
to appear, neglecting the flow of the marginal couplings, can 
now be obtained with the help of refermionization procedure 
as described in Ref. |l_3|. Specifically, the model ([8]) with two 
competing relevant couplings g\ and 174 maps onto the theory 
of four (weakly interacting) Majorana fermions which are or- 
ganized into a triplet with mass m t oc 2irv(Gi — G4) and a 
singlet with mass m s oc — 2irv(3Gi + G4). The Haldane-to- 
CD transition corresponds to the vanishing of the triplet mass, 
m t = 0, and takes place when J± — 2J X + 5J 2 /(tt 2 J) = 0. 
The CD-to-RS transition is described by the condition m s = 
and results in J± — 2J X + J 2 I {^J) = 0. Putting these two 
boundaries together leads to the estimate (|5). 

The neglect of the marginal terms, and in particular of the 
inter-chain current-current one, G 2 , is well justified in the 
asymptotic limit Ji. x / J — > 0. Away from this limit, but still 
keeping Jx t %/ J 1, one has to be mindful of the marginally 
relevant character of G 2 term in the case of antiferromagnetic 
(positive) sign of exchanges J± yX - Left alone, G 2 would di- 
verge at£ ~ J/J x , as follows from (fT3l . It may happen that 



for not too small J± yX /J the marginal coupling G 2 reaches 
value of the order 1 before the relevant couplings Gi,4 do so. 
In the following we denote the corresponding RG scale as £3: 
thus G 2 (4) 1 while G M (4) < 1. 

Such a behavior implies direct first-order transition between 
Haldane and RS phase a 13 i 19 ~ — , which can be analyzed with 
the help of abelian bosonization along the lines sketched pre- 
viously in Ref 1 1 3l In terms of conjugate bosonic fields <p m and 
6 m , where m = ± denotes symmetric/antisymmetric combi- 
nation of the original chain fields, the interchain interaction 
([8]l acquires the form 




2gi cos[\/47r#-] — (<7i — 34) cos[V47r^ + ] (16) 
+ {91 + .94) cos[\/47n ( 2_]^ . 

This expression includes only relevant and marginally rele- 
vant terms of ([8]) (i.e. terms N, E and J) evaluated at £ — £3. 
The two remaining terms (M and K) are omitted because of 
their irrelevancy at this late stage of RG development. In the 
regime we are interested here g 2 ^> 31 4 which implies that 
H! Dtei is minimized by (a) co + = y / 7r/2,0_ = 0, and (b) 
ip+ = 0,9- = y/Tr/2. The first choice describes Haldane 
state while the second corresponds to the RS one, see Reffl3l 
Since is pinned in both states, the last term in ( fTBI l ef- 
fectively averages to zero, cos[\/47r<yS_] — > 0. Treating the 
remaining two terms as a perturbation, we find that 34 = 3<?i 
describes a first-order transition line separating Haldane and 
RS phases. Everywhere on this line the energies of the two 
states are equal. Line's endpoints can be estimated as points 
where g 2 = 31,4, see Ref 1 1 31 and Figure 1 therein. Since the 
Haldane and the RS states are degenerate on the first-order 
transition line, elementary excitations on this line are (mass- 
less) domain walls or kinks, interpolating between vacua of 
types (a) and (b) above. The spin of such a kink can be evalu- 
ated as 

S t z otal = SI + 5| = d x (^ + ^2)1^ = d x ip+/V^- (17) 

Assuming, for example, (p+(x = —00) = and <p+(x — 
+00) = y/ir/2, we find 5 t z otal = 1/2: the kink is a spinonl 21 
These spinons can be visualized as spin- 1/2 end states of the 
Haldane phase segments - since the length of the segment is 
arbitrary, the spinons are mobile and massless 13 . Away from 
the line 34 = 3gi the spinons are bound into pairs as the en- 
ergy cost of two distant kinks is proportional to the length of 
'minority' phase segment between them. 

Once the initial parameters are such that one of the rele- 
vant couplings 314 reaches 1 before the marginal <? 2 , the in- 
termediate phase separating Haldane and RS phases becomes 
unavoidable. For positive 34 such intermediate phase real- 
izes staggered dimer (SD) phase (</?+ = \/tt/2, if- = y/ir/2) 
while negative 34 results in the CD phase (ip + = 0, ip- = 0) 
which is the focus of this work. Such a situation is realized 
in the strict weak-coupling limit, J±. x /J —> 0, when the 
relevant couplings are certain to flow to values of order one 
well before any of the marginal ones can do so - and thus we 
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once again conclude that dimerized phase is unavoidable. Its 
numerical detection is however very challenging in this very 
limit. 

Solving full system of RG equations (fT2l) numerically we 
find that conditions favoring the direct first order transition are 
realized for (approximately) J x > JJ s» 0.3 and J± rj 2J x . 
Even though it is not a priori clear that our weak-coupling the- 
ory can adequately describe these somewhat large inter-chain 
couplings, observation of the direct transition between Hal- 
dane and RS states in several numerical studies suggests that 
analytical description based on dT2l > and ( fToT l remains valid 
there. We thus conclude that antiferromagnetically coupled 
frustrated ladder should exhibit CD phase in the range ap- 
proximately given by <(3j and for not too large inter-chain 
exchanges, J x < J* . Stronger inter-chain exchange leads 
to the collapse of the CD phase and direct first-order transi- 
tion between Haldane and RS ground states. The situation is 
very similar to the phase diagram of the extended Hubbard 
model, 22 24 where one finds that charge-density-wave and 
spin-density-wave phases are separated by the bond-charge- 
density-wave phase. This intermediate phase collapses and 
gets replaced by the line of direct first-order transition once 
Hubbard repulsion constant exceeds some critical value. 

The outlined reasoning suggests simple way to avoid the 
marginal interaction issue altogether: all we need to do is to 
change the signs of both inter-chain exchanges to the ferro- 
magnetic (negative) ones. This simple change preserves frus- 
trated nature of the inter-chain couplings - large negative J± 
leads to an effective spin-1 chain and the Haldane phase while 
large negative J x forces rung spins into the RS phase. Impor- 
tantly, this change makes inter-chain G2 marginally irrelevant 
so that G 2 (l) = G 2 (0)/(1 + \G 2 (0)\£) -> -1/1 (neglecting 
for a moment effect of all other couplings). In effect, this sim- 
ple sign change sends the scale £3 to infinity, £3 — > 00. We 
now should be able, within the weak-coupling approximation, 
to get rid of G2 which forces the CD phase to shrink, without 
suppressing the all-important competition between G\ and G4 
terms. This is indeed observed in numerical solution of the 
system ( TTSI i. The difference between antiferromagnetic and 
ferromagnetic inter-chain couplings is illustrated in Figs. Q] 
and|2] Coupling G2, shown by dotted (green) line there, over- 
takes relevant G1.4 on the way to strong coupling in Fig. Q] 
(antiferromagnetic case) while it remains small in Fig. ^fer- 
romagnetic case). 

Results of numerical solutions for some sample values of 
J x . j_ are summarized in the TableU which shows that making 
the inter-chain exchanges ferromagnetic indeed helps one to 
unmask the novel columnar dimer phase. 
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FIG. 1: (Color online) RG flow of couplings 01,2,4,5 for the case of 
antiferromagnetic inter-chain exchanges J x = 0.15, J± = 0.296. 
Notations are as follows: Gi (red/dashed), G2 (green/dotted), — G4 
(blue/solid), Gg (magenta/dot-dashed), t denotes RG scale. 
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FIG. 2: (Color online) RG flow of couplings Gi.,2,4,5 for the case of 
ferromagnetic inter-chain exchanges J x = — 0.15, Jx = —0.306. 
Notations are as follows: —G\ (red/dashed), — G2 (green/dotted), 
— G4 (blue/solid), G5 (magenta/dot-dashed). £ denotes RG scale. 



results, we determine the ground-state phase diagram, which 
includes the CD phase in a wide parameter range between the 
Haldane and RS phases. 



III. FERROMAGNETIC INTER CHAIN COUPLINGS 

Motivated by the result of RG analysis in Sec. [II] we study 
the case of ferromagnetic inter-chain couplings J± , J x < 0. 
First, we treat the limit of strong rung coupling, |Jj_| 3> 
J, I J x |, and show that for J + J x =0 the model exhibits the 
CD long-range order. We then present our numerical DMRG 
and exact-diagonalization data for the model. Combining the 



A. Strong rung-coupling limit 

We consider the limit of strong rung coupling, |Jj_| ^S> 
J, I J x I • We first diagonalize the rung Hamiltonian H lung , 
whose ground states are a direct product of triplet states in 
each rung. We then include the effect of iJi cg and i/diag per- 
turbatively. It is convenient to rewrite the perturbation term 
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TABLE I: Brief summary of numerical solution of RG system d!2t 
for several values of J x . Positive (negative) values of J x ,x cor- 
respond to antiferromagnetic (ferromagnetic) inter-chain exchanges. 
Second (third) column describes the range of J± for which the 
columnar dimer phase is realized according to Eq.J 12t (Eq.lO). Last 
column shows the most divergent coupling constant which reaches 
value of order one first. 



Jx 


range of J± (eq.dl2t) 


estimate of J± (eq.(|5j) 


leading g 


0.05 


(0.09945, 0.09955) 


(0.0987, 0.0997) 




0.1 


(0.1978, 0.1982) 


(0.1949, 0.1990) 


.94 


0.15 


(0.295, 0.2956) 


(0.2886, 0.2977) 


.94 


0.2 


(0.39096, 0.392) 


(0.3797, 0.3959) 


.94 


0.3 


none 


(0.5544, 0.5909) 


.92 


0.4 


none 


(0.7189, 0.7838) 


.92 


0.5 


none 


(0.8733, 0.9747) 


92 


-0.15 


(-0.303, -0.3062) 


(-0.3023,-0.3114) 


94 


-0.3 


(-0.608, -0.635) 


(- 0.6091, -0.6456) 


.94 


-0.4 


(-0.815, -0.882) 


(-0.8162, -0.8811) 


.94 


-0.5 


(-1.02,-1.16) 


(-1.025, -1.127) 


94. 



as 

H = H\ c 
1 



H, 



(J + J x ) / v (Sl, n + ^2,n) • (Sl jn +1 + S2,n+l) 

n 

+ ^(^ ~ J~x) ^~^(Sl,n — S2,n) ' (Sl,n+1 — S2,n+l)- 

n 

(18) 

Note that the first term preserves the total spin in each rung 
while the second term changes the rung-triplet state to rung- 
singlet one and vice versa. 

When J + J x 7^ 0, the calculation is easy. The first term 
in Eq. ( fT8l gives a nonzero contribution at the first order per- 
turbation and lifts the ground state degeneracy of H lung . The 
effective Hamiltonian turns out to be the spin-1 Heisenberg 
chain, 



(19) 



where S n is the spin-1 operator consisting of rung spins Si,„ 
and S 2 ,„ and jW — (J + J x )/2. Therefore, if J + J x > 0, 
the system is in the Haldane phase, while the system exhibits 
the ferromagnetic ground state for J + J x < 0. 

For J + J x =0, the first-order perturbation vanishes, and 
we must turn to the second order. From a straightforward 
calculation, we obtain the second-order perturbation Hamil- 
tonian of the form, 



#(2) = J(2) 



- 1 



with 



J(2) 



"(J- Jx 



21 J 1 



2 Ji 



(20) 



(21) 



Therefore, the low-energy physics of the system is described 
by the spin-1 pure biquadratic chain with negative j' 2 '. For 
this case it has been established that the model has the dimer- 
ized ground state^^ii Hence, mapping the spin-1 dimerized 
phase back to our model, we conclude that the spin- 1/2 two- 
leg frustrated ladder (HJ must exhibit CD phase along the line 
J x = — J in the strong rung-exchange limit. 



B. DMRG results 

To search for the CD state and determine the ground-state 
phase diagram, we carry out the DMRG calculation 3 ^ for 
the frustrated ladder ([T). The calculation is performed for the 
system with up to L = 192 rungs. For the efficiency of the 
DMRG method, the open boundary condition is imposed in 
the calculation. The number of kept states are typically m — 
350 for L < 96, m = 400 for L = 128, 192, and up to 
m = 500 for some cases of the severe truncation error. We 
have monitored the truncation error of the data by comparing 
the results obtained with different rn's and confirmed that the 
m convergence has been achieved for the data shown in the 
following. 

To detect the CD order, we calculate the local CD operator 
in the open ladder with L rungs, 

DcT>{n;L) = ((Sj, n • Sj :7l+ i) ~ (Sj in+ i ■ S^ n+2 )) , 
3=1,2 

(22) 

where (• • • ) denotes the expectation value in the ground state, 
i.e., the lowest-energy state in the subspace of zero magne- 
tization 5 t z ot = • S z j n = 0. In the CD phase, the CD 
order induced at open boundaries of the ladder penetrates into 
the bulk and exhibits a long-range order. In the other phases 
with a spin gap, the CD order is expected to decay exponen- 
tially when we move from the boundary into the bulk, while 
we expect that the CD order decays algebraically at a critical 
point. We may therefore be able to identify the CD phase and 
transition points by monitoring the system-size dependence 
of the CD operator Dcv(n; L) at the center of open ladder, 
n = L/2. In the calculation, we set L to be a multiple of four 
so that Dqd{L/2\ L) is positive. 

Figure [3] shows the dependence of the CD operator 
Dcb(L/2; L) on the rung coupling J± for J = 1 and sev- 
eral fixed J x . We find that Dcb(L/2; L) has a broad peak, 
indicating that the CD order is strong in a rather wide regime 
of J±. We note that for J x = — 1 and J± < —2, the L con- 
vergence of Dcd(£/2; L) seems almost achieved, suggesting 
the appearance of the CD long-range order. 

In order to determine whether or not the CD order sur- 
vives in the thermodynamic limit, we investigate the system- 
size dependence of the CD operator Dcb(L/2; L). Figure [4] 
shows the L dependence of the CD operator Dcn(L/2; L) 
for some typical sets of coupling parameters. It is clear that 
Dcb(L/2; L) for J x = —1 and large negative J± converges 
to a finite value at L — >• 00. 

We also show in Fig. [5] the spatial profile of the nearest- 
neighbor spin correlations (Sj >n ■ Sj. n +i) — C av + (— l) n (e n ) 




—1 








CM 


0.4 












0.6 


Q 






0.4 


O 






Q 


0.2 


_ -j 








Q 
O 


0.2 






Q 





(c) 



FIG. 3: (Color online) CD operator at the center of the open ladder, 
Dcd(L/2; L) as a function of Ji for J = 1 and (a) J x = —0.2, 
(b) J x = —0.5, and (c) J x = —1.0. The symbols represent the 
data for L = 16, 24, 32, 48, 64, 96, 128 and 192 from top to bottom. 
Inset in (c) shows the data for J x = —1.0 and broader regime of J±, 
-8 < J± < 0. 
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for (J, J x , Jx) = (1, —1, —3), which clearly demonstrates 
the presence of well-developed columnar dimer order. The 
average energy density, calculated in the middle of the lad- 
der, is found to be C av = ((S jVL/2 _ 1 ■ S jVL/2 ) + (S jVL/2 • 
£>j,l/2+i))/2 = —0.384. We see strong modulation of the 
bond energy (Sj tTi ■ Sj, n +i) between even and odd bonds. 
The amplitude of the modulation saturates in the middle of 
the ladder where the bulk dimerization value is achieved, 
(e n ) — > 0.136. We find that bond modulations in the two 
chains are in-phase, implying columnar ordering of stronger 
and weaker bonds. This finding represents direct proof of the 
CD phase in the frustrated ladder model (fl]i with ferromag- 



FIG. 4: (Color online) System-size dependence of the CD operator 
at the center of the open ladder, Dcd {L/2; L), in a log-log scale for 
J = 1 and (a) J x = -0.2, (b) J x = -0.5, and (c) J x = -1.0. 



netic inter-chain exchanges. 

For smaller | J x |, on the other hand, the appearance of the 
CD long-range order is not so clear; Dcr>(L/2; L) still de- 
creases with L even at the largest L calculated [see Fig. [4] 
(a) and (b)]. However, we find that in some parameter regime 
Dct>(L/2; L) bends upward in a log-log plot. This means that 
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FIG. 5: Spin correlations between nearest-neighbor spins (bond en- 



ergy) along the legs of the ladder, (Sj !n ■ Sj,„+i), for (J, J x 



(1, -1, -3) and L ■. 
j = 2 are identical. 



192. The correlations in the legs j 



, Jx) = 
= 1 and 



the decay of Dcv{L/2; L) becomes slower as L gets larger, 
which suggests the emergence of the CD long-range order in 
the thermodynamic limit. 

To elucidate the bending-up behavior, we also investigate 
the system-size dependence of the slope of the log-log plot, 



A CD (^i) 
\og[D C D{Li 



-i/2;Li + i)]-log[DcD(£i/2;£< 



log(it 



log(ii) 



(23) 



where Xi — (Li + Li + i)/2 and Li = 
16,24,32,48,64,96,128,192 for i = 1,2,..., 8. If 
Dcd(L/2; L) decays exponentially with increasing L, the 
slope Acd(^) decreases as x increases. If Dcd(L/2; L) 
exhibits a long-range order, Acd(^) increases with x and 
converges to zero at x oo. Furthermore, if Dcd(L/2; L) 
decays algebraically, Acd(^) converges to a finite negative 
value at x — > oo. Figure [6] shows the data of Acd(^) as 
a function of The results clearly suggest that there are 
parameter regions where Acd(^c) increases with x. We take 
this behavior as an evidence of the CD phase. 

Based on the above results we conclude that the CD phase 
emerges in a finite region in the J± - J x plane. The phase 
boundaries estimated from the results of the slope Acu(x) 
above are plotted in the phase diagram, see Fig.|9]in Sec. lIIIDl 
We note that, as shown in the Appendix, the bending-up be- 
havior of the dimer operator in the log-log plots is also ob- 
served in the frustrated Heisenberg chain (which can also be 
viewed as the zigzag ladder), which is well known to exhibit 
the dimer phase for sufficiently large next-nearest-neighbor 
exchange J 2 ^— This observation provides us with an im- 
portant check of the approach to the frustrated ladder (Q~|i and 
supports our interpretation of the data in Figs.|4]and|6] 

Although not expected from the RG analysis, we have also 
examined possibility of the SD order in the model. For this 
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purpose, we have calculated the local SD operator, 

D S n(n;L) 

= ^2 {(Sj,n ■ Sj.n+l) - (Sj>+1 ■ S J>l+2 )) 



3 = 1,2 



(24) 
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FIG. 7: (Color online) System-size dependence of the SD operator 
at the center of the open ladder, Dso{L/2; L), in a log-log scale for 
J = 1 and (a) J x = -0.2, (b) J x = -0.5, and (c) J x = -1.0. 



in the frustrated ladder ([T]) with up to L = 96 rungs. In the 
calculation of the SD operator, we have employed an open 
boundary condition with an extra spin at each edge, which se- 
lects one of the SD patterns and lifts the two-fold degeneracy 
in the possible SD ground states. The results are presented 
in Fig. |7] We find that D^(L/2\ L) bends downward in a 
log-log plot, indicating the exponential decay. [Dgu(L/2; L) 
for J x = —0.2 and J± = —0.41, at which point we have 
found that the decay of the SD order is the slowest, exhibits 
a nearly-linear behavior, but it actually bends down slightly.] 
We have performed the same calculation for a wide parameter 
regime and found that Dsb(L/2; L) decays exponentially in 
each gapped phase or, at most, decays algebraically at a tran- 
sition point. We thus conclude that the SD phase is absent in 
the model (Q~|i with ferromagnetic J± and J x . 



C. z operator 



Here, we discuss another numerical approach to the prob- 
lem, based on so-called "z operators" which are used to 
distinguish different valence-bond-solid (VBS) states in one- 
dimensional spin systems. For the frustrated ladder model ([T), 



^diag(i) 



(exp 



(exp 



2tt 
2tt 



n=l 
L 



(25) 



It has been shown 39,40 that the z operators in the spin- 1/2 two- 
leg ladder with L rungs under the periodic boundary condition 
exhibits the following asymptotic behavior with L, 



^■ung/diagW ~ (-l)^™ 3 [1 - 0(1/ L)] 



(26) 



where AVbs is an integer depending on the VBS pattern of 
the state under consideration: it represents the number of sin- 
glet bonds 'cut' by a line parallel to the rung/diagonal link. 
The z operator then measures topological parity of the dimer 
covering pattern describing particular gapped state. In our 
case, z rung converges to 1 for the RS and CD states (even 
number of singlets crossed) in the thermodynamic limit, while 
Zrung — > — 1 for the Haldane state (the number of crossed sin- 
glets is always odd). Conversely, Zd 



lag 



-1 for the RS and 
CD states, while Zdiag — > 1 in the Haldane state. A remark- 
able feature of the z operators is that they change their sign at 
the transition between phases having different parity of A^vbs • 
This property makes the z operators more powerful in deter- 
mining the critical point of such a phase transition than the 
string order parameter, which just vanishes at the transition. 20 
Indeed, the z operators have turned out to be successful in de- 
termining the direct RS-Haldane transition point occurring for 
large antiferromagnetic Jj_ iX F° r tne present case of ferro- 
magnetic Ji jX , we can use z rung and Zdi ag to locate the tran- 
sition point between the CD and Haldane phases. 

Using the exact-diagonalization method, we have calcu- 
lated the z operators, z rung and Zdiag, in the ladder (HJ with 
up to L = 12 rungs under the periodic boundary condition. 
Figure [8] presents the results for typical parameter lines with 
J = 1 and fixed J x . For J x < 1, we have observed the 
sign change in z r ung (Zdiag) from positive (negative) to neg- 
ative (positive) values as J± decreases. The crossing point 
of z rung and Zdi ag thereby gives an estimate of the transi- 
tion point between the CD and Haldane phases. While the 
L dependence is negligibly small for small | J x |, the crossing 
point for large | J x | moves sizably with L, suggesting that the 
finite-size effects still remain. However, we emphasize that 
the crossing point shifts towards smaller J± with increasing 
L, which means that the range of CD phase broadens as L 
increases, and approaches smoothly to the CD-Haldane tran- 
sition point obtained from the DMRG analysis. (See also the 
phase diagram, Fig.|9]in Sec. lIIlDl ) Thus, we can safely state 
that the analysis of z operators also supports the appearance of 
the CD phase. For J x = 1, z rung (zdiag) is positive (negative) 
for the entire regime of J± calculated. The result is consistent 
with the prediction of the perturbative analysis in Sec. IIII Al as 
well as the DMRG results in Sec. IIII Bt which show that the 
CD phase extends to the limit J± —> — oo. 
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FIG. 8: (Color online) J± dependence of the z operators for J = 1 
and (a) J x = -0.2, (b) J x = -0.5, (c) J x = -0.8, (d) J x = 
— 1.0. Dark (blue) and light (red) curves represent z IU n S and 2dia g , 
respectively. 



D. Phase diagram 



FIG. 9: (Color online) Ground-state phase diagram for J — 1 and 
the ferromagnetic interchain coupling, J± < and J x < 0. Open 
circles represent the transition points obtained from the DMRG cal- 
culation while the other symbols show the CD-Haldane transition 
points from the analysis of z operators. Dotted lines are guide to 
eye. In the schematic pictures for each phase, the bold lines repre- 
sent singlet pairs and the ellipses stand for the symmetrization of two 
spins. 



Combining the above results, we determine the ground- 
state phase diagram in the parameter plane for ferromagnetic 
J_l and J x . Figure [9] shows the resultant phase diagram, 
which includes the Haldane, RS, and CD phases. We clearly 
see that the CD phase appears in a wide parameter region, 
which is seen to expand as \J±\ and | J x | become bigger. 
The transition line between the RS and CD phases seems to 
nearly coincide with the line of J± = 2 J x . The boundary be- 
tween the Haldane and CD phases starts from J± = J x = 
and runs towards smaller J± as J x decreases, approaching 
smoothly the limit of the strong rung-exchange, J x = — J 
at J± — > — oo. It is worth noting that the DMRG result on 
the RS-CD transition line agrees even quantitatively with the 
result of RG analysis in Table U and the behavior of the CD- 
Haldane transition line is also consistent with the analytical 
RG result. This observation strongly supports the correctness 
of the RG analysis in SecHIl 



IV. ANTIFERROMAGNETIC INTER-CHAIN COUPLINGS 

The numerical results in Sec.|lII|have revealed that the frus- 
trated ladder ((TJ with ferromagnetic J± and J x exhibits the 
CD phase in a wide parameter regime, in agreement with the 
prediction of RG analysis in Sec.|ll] Since the validity of the 
RG analysis relies only on the small amplitudes of the inter- 
chain couplings J± and J x and is not affected by their signs, 
we naturally expect that the RG analysis is correct also for 



the antiferromagnetic couplings. To examine the expectation, 
we re-visit the frustrated ladder (fl} with antiferromagnetic J± 
and J x . For this case, it has been shown rather clearly that for 
large J± and J x the direct first-order transition takes place 
between the RS and Haldane phases, 15 while the situation is 
still controversial for small J± and J x To clarify the sit- 
uation we have performed the DMRG calculation for a param- 
eter line J = 1 and J x = 0.2 and investigated behaviors of 
the CD and SD operators. 

Figures [10] and QT| show the system-size dependence of 
the CD and SD operators at the center of the open ladder, 
Dct>{L/2\ L) and Dgr>(L/2; L), and the slopes of their log- 
log plots, Acd(^) and Asd(x), respectively^!! [Asd(x) 
is defined in the same way as Eq. d23ll.l We note that our 
data of the CD operator D C d(L/2;L) for J± < 0.37 and 
J± > 0.39 coincide with the results shown in Ref. 14, while 
the data for J± = 0.38 was not presented there. We find 
in Fig. [TO] that both D C d{L/2;L) and D SD (L/2;L) decay 
exponentially with L for J± < 0.37 (Haldane phase) and 
J± > 0.39 (RS phase), suggesting the absence of the CD and 
SD orders in the parameter regions. On the other hand, it is 
remarkable that for J± = 0.38 the CD operator Dcd(L/2; L) 
bends upward in the log-log plot, indicating the emergence of 
the CD long-range order. The tendency toward the CD order- 
ing is elucidated also in Fig. II U a). which shows that the slope 
of the log Dco(L/2; L)-\ogL plot, Acd(x), increases with 
x. We note that, in contrast to the CD operator, the SD oper- 
ator Dsd(L/2; L) exhibits the bending-down behavior in the 



10 



Q 

O 

Q 



10" 







Jx= 0.2 


-o- 


Jx= 0.35 \ 






Jx= 0.37 \ 




-v- 


Jx= 0.38 




-a- 


Jx= 0.39 




-<- 


Jx= 0.41 





in agreement with recent numerical finding in Ref.[l6l 



Iff 5 



Q 



10" 







Jx= 0.2 


-o- 


Jl= 0.35 ^\ 






Jx= 0.37 \ 






Jj_= 0.38 




-a- 


Jx= 0.39 






Jx= 0.41 





(b) 



iff 5 



FIG. 10: (Color online) System-size dependence of the dimer op- 
erators at the center of the open ladder in a log-log scale for J = 1 
and J x = 0.2; (a) the CD operator D C o(L/2; L) and (b) the SD 
operator Dsd{L/2; L). 
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FIG. 1 1 : (Color online) Slope of the log-log plots of the dimer oper- 
ators for J = 1 and J x = 0.2; (a) Acd(^) and (b) Asd(s). 



log-log plot even for Jj_ = 0.38. The opposite trends of the 
CD and SD operators imply that the growth of the CD order 
observed at J± = 0.38 is not a critical enhancement at a tran- 
sition point but an indication of a true CD long-range order. 
We therefore expect that the CD phase appears in a narrow 
but finite parameter region around J± = 0.38, in accordance 
with the RG prediction 13 and the discussion in Section|II] and 



V. DISCUSSION 



The main result of our study is the discovery of the colum- 
nar dimer phase in the frustrated ladder problem with fer- 
romagnetic inter-chain interactions, see Fig. [9] This find- 
ing, confirmed by extensive DMRG analysis in Sec. IIII Bl 
is based on analytic RG arguments summarized in Sec. [TT] 
It confirms novel mechanism of dimerization by frustrated 
interchain couplings, proposed in Ref. 13. Previous sight- 
ings of the spontaneously dimerized state, of either columnar 
or staggered type, were restricted to models with four-spin 
interaction terms, such as the ring-exchange model and the 
SU(2)xSU(2) ladder.^ii 

The success of this study in describing ferromagnetic inter- 
leg exchanges gives us confidence in essential validity of the 
weak-coupling RG approach and makes it possible to re-visit 
the more complicated case of antife rromagnetic inter-leg ex- 
changes, as described in Sec. HV] There we also find hints 
of developed CD order at (J, J x , J x ) = (1,0.2,0.38), in 
agreement with Ref. 16. The extent of the CD region is very 
narrow: finite-size scaling analysis in Ref. [l6j estimates that 
0.373 < Jj_ < 0.386 for J x = 0.2. Such a limited range may 
explain negative results of the two previous studie s - 14 ' 15 . 

In addition to these numerical observations our work takes 
important step forward in uncovering the reason for the more 
narrow than naively expected, on the basis of the estimate ©, 
range of existence of the CD order. That feature, as we ar- 
gue in Section|IIl has to do with marginally relevant character 
of the current-current interaction between spin chains in the 
case of antiferromagnetic inter-leg exchanges. We predict that 
the CD phase ceases to exist at all once inter-leg exchange 
J x exceeds the critical value of the order 0.3. Connecting 
this CD phase with the dimerized phases of frustrated two- 
dimensional spin models (see Ref. |48| for the original large-N 
study and Ref. |49| for recent developments) represents an im- 
portant outstanding problem. 

Before concluding we would like to note that there exists 
another simple route to the dimerized phase. It consists in 
turning marginally irrelevant in-chain backscattering G5 into 
a marginally relevant one. 50 This is achieved by introduc- 
ing sufficiently strong antiferromagnetic coupling J2 between 
next-nearest spins along the legs of the ladder. Provided that 
it exceeds the critical value 17 , J2 > 0.241J, the legs of the 
ladder will be spontaneously dimerized even in the absence of 
any inter-chain coupling. The remaining weak inter-chain in- 
teractions then work to stabilize one of the two ordered dimer- 
ization patterns, columnar or staggered, as is described in Ref. 
50 and observed in Ref. IB Connecting this large- J2 regime 
with the case studied here represents another interesting topic 
we leave for future. 
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Appendix: Dimer order in J1-J2 chain 

In this Appendix, we check the behavior of the dimer opera- 
tor in a finite open spin chain as a function of chain length. To 
this end, we consider well-understood frustrated Heisenberg 
chain (Jx - J2 model), 

H Z i g = J\ S„ ■ S„ + i + J2 S n • S„ + 2, (A.l) 

n n 

where S„ is the spin- 1/2 operator at the nth site and J\ 
and J2 are coupling constants of the nearest- and next- 
nearest-neighbor exchange interactions, respectively. It is 
well established that in the case of antiferromagnetic cou- 
plings, Jx, J 2 > 0, the J\-Ji chain dA.lt exhibits a criti- 
cal (Luttinger-liquid) phase for J2/J1 < (Ja/Ji)c — 0.241, 



while for J%l J\ > (J2 /Ji ) c the ground state is spontaneously 
dimerized.-ii^— 

Using the DMRG method, we calculate the dimer operator, 

D(n; L) = (S„ • S n+ i - S„+i • S n+2 >, (A.2) 

in the chain with up to L = 128 spins under the open bound- 
ary condition. Figure [12] shows the system-size dependence 
of the dimer operator at the center of the chain, D(L/2; L), 
and its slope in the log-log plot, A(x), for several typical val- 
ues of J2/J1. [The slope A(x) is defined as in Eq. d23l ).l For 
J2/J1 < 0.241, where the model is in the critical phase, the 
dimer operator D(L/2; L) decays algebraically with L, as ex- 
pected. For 0.241 < J2/J1 < 0.3, for which regime it is 
known that the system is in the dimer phase but the spin gap 
is exponentially small, the dimer operator seemingly decays 
in a power law. This can be understood as a consequence of 
the fact that the correlation length is so large that we can not 
reach the asymptotic behavior of the dimer operator within 
the system size treated, L < 128. For J2/J1 > 0.35, deep 
in the dimer phase, D(L/2; L) shows the bending-up behav- 
ior in the log-log scale, and eventually, the dimer long-range 
order is clearly observed for J2/J1 = 0.45. 

The results indicate that the bending-up behavior of the 
dimer operator in the log-log plot is observed only in the dimer 
phase and when the system size is comparable to or larger 
than the correlation length. We can therefore safely regard the 
bending-up behavior as an evidence of the dimer ordering. 
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